Efficiency of autonomous soft nano-machines at maximum power 
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We consider nano-sized artificial or biological machines working in steady state enforced by im- 
posing non-equilibrium concentrations of solutes or by applying external forces, torques or electric 
fields. For unicyclic and strongly coupled multicyclic machines, efficiency at maximum power is 
not bounded by the linear response value 1/2. For strong driving, it can even approach the ther- 
modynamic limit 1. Quite generally, such machines fall in three different classes characterized, 
respectively, as "strong and efficient", "strong and inefficient", and "balanced". For weakly cou- 
pled multicyclic machines, efficiency at maximum power has lost any universality even in the linear 
response regime. 
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Introduction.- Molecular motors [H as well as the 
recently developed artificial nano-machincs inspired by 
them 0, Q operate in an aqueous solution of constant 
temperature. In contrast to heat engines limited by 
Carnot's law, thermodynamics constrains their efficiency 
by 1. Like for heat engines, however, operating at the 
upper bound comes at the price of zero power since 
it requires infinitely slow driving. A practically more 
relevant question then is about efficiency at maximum 
power (EMP). For heat engines this problem has a few 
decades of history [1] where the full picture exhibiting 
some universality close to the linear response regime has 
just emerged [5j— lsf| . 

For machines working under isothermal conditions, 
which provide arguably a much more realistic route for 
a successful implementation on the nano to micro scale 
than dealing with heat baths of different temperature, 
efficiency has been discussed mostly in the context of 
molecular motors (ami. The issue of EMP, however, has 
received much less attention so far. For molecular mo- 
tors, a case study has shown values well above the linear 
response result 1/2 17(. On the other hand, this value 
has recently been derived as a universal bound even be- 
yond the linear response regime under conditions claimed 
to be relevant for mesoscopic machines with a sufficiently 
large number of internal states [l8| . 

In this paper, we address the issue of EMP for au- 
tonomous soft nano-machines using a quite general ap- 
proach requiring minimal assumptions. "Autonomous" 
stands for steady state conditions as typical for molec- 
ular motors and envisaged for artificial machines. At 
present, the latter mostly still need a modulation of ex- 
ternal parameters for driving them through a cycle as in- 
vestigated theoretically in Refs. [US!- "Soft" is short 
for working in an aqueous environment which will require 
some flexibility in the molecule (s) allowing for conforma- 
tional changes to operate as a machine. By implementing 
thermodynamic consistency from the very beginning and 
by separating universal from system specific quantities, 
the result will be applicable to a large class of machines 
which turn out to fall into three regimes. In particular, 



we find that EMP is not bounded by 1/2. It can rather 
approach 1, which, however, requires some care in the 
design of the machine as well as a sufficiently large pa- 
rameter space available for the maximization of power. 
In fact, for a severly restricted parameter space, EMP 
looses any universality. 

Model.- We model the kinetics of the nano-machine as 
a Markov process in a heat bath of constant tempera- 
ture 22I-424J] . At any time, the machine is in one of sev- 
eral possible states. A transition between state m and 
state n happens with a rate k mn . Typically, in a tran- 
sition, some (external) quantities d" m later to be associ- 
ated with the function of the machine change. Examples 
for such quantities could be (i) position along a linear 
track, (ii) rotation angle, (iii) number of consumed (or, if 
negative, produced) molecules of a certain species, or (iv) 
number of charges transported against an external field. 
For each such quantity, there exists a "conjugate" exter- 
nal field h a with the property that the product h a d^ in 
is a (free) energy. Specifically, for the above cases, the 
conjugate field is (i) force /, (ii) torque N, and (iii) de- 
viation A// 1 of chemical potential of species i from its 
equilibrium value and (iv) potential difference A(f>. 

For the object to operate as a useful machine, at least 
two of the fields, an input field h m and an output field 
h out , have to be set to non-zero values. For a rotary 
molecular motor like the Fl-ATPase (see, e.g. 25|, l26j 
and Refs. therein), one has to provide ATP at higher 
than equilibrium chemical potential, i.e., h ln = A/x ATP . 
Hydrolysis of ATP can then be used to move against an 
external torque, i.e., h = N. In contrast to macro- 
scopic machines, the role of input and output is not 
unique but can be interchanged depending on the in- 
tended "purpose" of the machine. The ATPasc can 
deliver mechanical work if fed with an excess of ATP 
molecules but can also synthesize ATP if pulled by an 
external torque in the opposite direction. 

The external fields will affect the transition rates. 
Thermodynamic consistency requires a local detail bal- 
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ance (LDB) condition of the form 



exp(/i in C„ + h^dTn ) (1) 



where (k mn /k nm ) eq refers to the equilibrium ratio where 
all external fields are set to zero. Here, and throughout 
the paper, we measure energies in thermal units (ksT = 
1). For constant external fields, h ln and h° ut , the machine 
reaches a steady state, i.e., will operate autonomously. 

Unicyclic machines.- We first discuss unicyclic ma- 
chines, for which the m = 1,...,M possible states are 
aligned in one cycle such that each state has two neigh- 
boring states. If the machine steps through the com- 
plete cycle in forward direction, it has consumed the 
work or (free) energy w m = J2m=i ^'"Xn.m+i an< ^ de- 



livered the output w 
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,1 h out d^ m+1 where 
^MM+i = ^mi- Likewise, had it stepped through the 
complete cycle in backward direction, it would have re- 
leased —Wi n to the input reservoir and consumed — w ou t 
from the output reservoir. In order to work as a machine 
in the conceived sense, the mean time r + it takes to 
complete the cycle in forward direction has to be smaller 
than the mean time t~ for completing the cycle back- 
wards. Both times can be expressed diagrammatically in 
terms of the transition rates but the explicit expressions 
will not be needed here [HI]. The key point is that the 
ratio r + /r~ is given by the ratio between the product of 
all backward rates and the product of all forward rates. 
With the local detail balance condition, cq. (JT|), and the 
identifications of input and output work just given, this 
ratio becomes 



r+/ T - = e b 



(2) 



Hence the power (or, more generally, the rate of "yield" 
if the output are chemical products) delivered by the ma- 
chine in the steady state is 



■Pout = w out (l/T + -1/t ) = Wout [1 - e v 



}/r + - (3) 



Likewise, the power used by the motor becomes 

P in = w in (l/r+ - 1/t-) = w in [l - e w ^- w -]/T+. (4) 

These transparent expressions for the power where the 
specific characteristics of the motor enter through the 
single quantity r + constitute our first main result. In 
the regime < w ou t < w; n , the machine will work as 
intended. Its efficiency is simply given by 77 = -P ut/-Pin = 
Wout/wi n and is obviously bounded by thermodynamics 
through < t] < 1. For w ou t = w- m , the machine has 
optimal efficiency 77 = 1 but does not deliver any power 
since it then cycles as often in forward as in backward 
direction. 

The concept of EMP requires to identify the admissible 
variational parameters. Rather than starting with an 



arbitrary parameter set {A^}, the form of eq. (|3|) strongly 
suggest to focus first on the following three choices: (i) 
Wout, (ii) Win, and (iii) w ou t and w; n . For each case, the 
crucial question becomes how the forward cycle time r + 
depends on input and output. We will first assume that 

x ou t = — d In t + /dw ou t and Xi n = — d\nr + /dw- ln 

are constants, i.e., that r + depends mono-cxponcntially 
(with cither sign) on input and output work as illustrated 
below with a specific example and later discuss the gen- 
eral case. 

Maximization with respect to output w ou t- _ For given 
input Win, the condition dP ou t/dw ou t = leads to the 
implicit relation 
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1 + x out w* 



(5) 



for the optimal output w* ut at fixed input Wi n - Conse- 
quently, we obtain for EMP under these conditions 
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(6) 

Two limit cases can be discussed analytically. 

First, for small deviations from thermodynamic equi- 
librium, Wi n <C 1, expanding (O and ((6]) leads to w* ut « 
Win/2 and to an efficiency at maximum power given by 

77 *nt - 1/2 + (*out + l/2)w in /8 + 0(wf„). 

For x ou t > — 1/2, EMP increases if the machine oper- 
ates beyond the linear response regime thus beating the 
bound 1/2 found in Ref. [18( under different conditions. 
Second, for large Wi n 3> 1, two cases must be distin- 



guished. If lout > 0, EMP, 77* 



1 - [lll( 1 + 1 /x ou t )] / Win , 



approaches the thermodynamic upper bound 1 . If x ou t < 
0, however, w ou t — > l/|£ ut| for Wi n — > 00 which im- 
plies 77* ut « l/|:EoutWin|. l n this case, efficiency at optimal 
power approaches with increasing input since maximum 
power is reached for a finite w* ut . 

Maximization with respect to input w- m .- We now as- 
sume that the output w ou t is given and search for the 
optimal input w* n (w ou t) > w ou t maximizing the power. 
Two cases must be distinguished, (i) For x m > 0, the 
power grows unbounded as w in increases and hence EMP 
vanishes for all w out . (ii) For x- m < 0, the optimal input 
reads w* n = w ou t + hi(l + l/|xi„|) implying the EMP 
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ln(l + l/|x in |) 



For Wout < 1, the efficiency n* n ps w out /ln(l + 1/1x^1) 
differs from the linear response result derived above for 
output maximization since the optimal w* n is always a 
finite distance from w ou t and hence no longer within 
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FIG. 1: Efficiency at maximum power r/**(xi n , x out ). Inset: 
Example of a multicyclic machine with 6 states, 6 cycles and 
the input transition (56). If (12) is the output transition (or, 
equivalently, (56) or (12)) the machine is strongly coupled 
(SC). For (23) as output transition, it is weakly coupled (WC). 
In the latter case, C\ consists of cycles (123456) and (12356), 
C 2 = {(2345), (235)}, C 3 = {(1256)}, and C 4 = {(345)}. 



the linear response regime. For io ou t 3> 1, one gets the 
asymptotics if m « 1 - [ln(l + \j\x- m \)\w oxlt . 

Maximization with respect to input and output.- This 
case, solved by w* ut and w**, leads to three regimes for 
EMP, n** ee w**/w** t , shown in Fig. ffl 

(i) For x ou t > max{— a;i n , 0}, the power at optimal out- 
put increases exponentially with growing input. In this 
case, one has w?* — > oo and w* ut — > oo with the finite 
difference wf* — w** t « ln(l + l/x out ) and hence n** = 1. 
Such machines could be called "strong and efficient" . 

(ii) For a; out < and x m > 0, the power still increases 
with increasing input, and, hence w n * — > oo. The op- 
timal output, however, in this limit remains finite at 
w out = Vl^outl which implies 77** = for these "strong 
but inefficient" machines. 

(iii) Finally, for x- nl < niin{— x outl 0}, the power peaks 
at a genuine maximum for w** t = l/[— (x- m + x ou t)] and 
*°*n = w out + ln(l + l/|a;i„|) leading to 
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t/ W in =!/[!+ I^out + Xin\ ln(l + l/\x h 



Such machines could be called "balanced" . 

Example. - For a simple but still instructive specific ex- 
ample, we consider a one-state molecular motor for which 
hydrolysis of ATP leads to a forward step of length d 
along a linear track in the direction of an applied force 
/ with rate k + = l/r+ = fc exp[A/i ATP + f6 + d], where 
ko is the equilibrium rate. A backward step with rate 



and k = 1/r = fc exp[A/i 



ADP 



A/r ! - f0-d}, ii 



volves synthesis of ATP from ADP and Pi. The load 
sharing factors 9 + and 6~ with 6 + + 8~ = 1 guarantee- 
ing the LDB condition eq. ((T|) are related to the distance 
of the activation barrier in forward and backward direc- 
tion, respectively [l3j |. The molecular architecture thus 



determines the sign of x ou t = d Wout lnk + = — 6 + since 
Wout = -fd. 

The input w in = A/i ATP - A^ ADP - Apt p ' can be 
changed in a variety of ways. The most obvious one, 
increasing A^ ATP , leads to X[ n = 1, whereas decreasing 
A/x ADP leads to a;; n = 0. If both chemical potentials are 
changed simultaneously, any value for be reached, 

in principle. Thus, depending on the molecular architec- 
ture, in the EMP diagram, Fig. [H this motor can either 
cross from balanced to strong and efficient (for 6 + < 0) 
or from balanced to strong and inefficient (for 9 + > 0) 
as the way the input power is delivered is changed. 

Arbitrary unicyclic machines.- So far, the complete 
analysis benefitted from using two simplifying assump- 
tions which will not always apply. First, we assumed the 
dependence of the forward cycle time r + to be mono- 
exponential in both input and output. While each indi- 
vidual rate will typically depend mono-cxponcntially on 
w out and Wi n , especially for intermediate values of w ou t 
and Wi n several rates may contribute comparably to r + . 
Then £ ou t becomes a function x ou t (w ou t , wi n ) where, e.g., 
for maximization with respect to output, w ut has to be 
determined self-consistently from eq. ([5]). 

Correspondingly, in the analysis of the limiting values, 
the quantities x ou t an d x- ln have then to replaced by their 
respective asymptotic limits. Gcncrically, in a unicyclic 
machine consisting of several states, increasing the input 
will speed up one (or several) forward transitions. If at 
least one forward transition is not affected by the increas- 
ing input, it will act as a bottleneck for the whole cycle 
and, hence, X[ n ~ for large w m . On the other hand, in- 
creasing the load, i.e., increasing w ou t will typically slow 
down at least one forward transition now becoming the 
bottleneck and consequently x ou t < 0. Thus, generically, 
according to the EMP diagram, Fig. [IJ such machines at 
maximum power will approach zero efficiency. In prac- 
tice, increasing w m will not lead to a significantly larger 
power since r + will approach a limit value for large W{ n . 
At any finite W{ n , one may thus still reach a reasonable 
efficiency as given by eqs. ([5]) and ([6]) and miss maximum 
power by only a small amount. 

The second major restriction so far was to focus on 
w- m and W7 ut as control parameters. For full generality, 
one should consider a set of arbitrary parameters If 
the admissible range of these parameters spans the whole 
sector < w ou t < w- m < °o (which will require at least 
two control parameters) then we are back at the case 
just discussed. On the other hand, if one has only one 
control parameter, or if only a limited part of the above 
sector can be coverered, EMP will become a strongly non- 
universal concept and its specific value will depend on the 
functional dependence of w ou t, wi n and r + on {A,}. 

Multicyclic machines.- In general, a nano-machine will 
contain several cycles. We will assume that the input 
affects only one transition, the "input transition" and 
that there is only one "output transition" which may 
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or may not be the same as the input transition. Two 
generic cases, strong coupling (SC) and weak coupling 
(WC) should then be distinguished, see inset of Fig. [1] 

(SC): If the input and the output transition arc iden- 
tical, or if there exists a direct connection between the 
two without any intermediate bifurcation, any cycle con- 
taining the input transition will also contain the output 
transition. For such strongly coupled machines exactly 
the same formalism as for unicyclic machines applies with 
the only caveat that r + appearing there is now given 



by l/i 



J2i^/ T i where the sum runs over all cy- 



cles that include input and output transition and the r 4 + 
are the corresponding forward cycle times. Thus, such 
strongly coupled machines obey the same relations for 
efficiency and EMP as discussed above for unicyclic ma- 
chines. Note, however, that as the explicit expressions 
would show the are affected by the rates contribut- 
ing to the cycles not containing the input and output 
transition. 

(WC): If the strong-coupling condition is not fulfilled, 
each cycle of the machine falls into one of four disjoint 
subsets C\, C2, C3, C4 containing, respectively, input and 
output transition, only the output transition, only the 
input transition, and neither one. Since only C\ and C2 
contribute to the output, the power is given by 



Pn 



«w[(l/r+ - 1/rf) + (1/t+ - l/r a -)] 



»-)/r+ - (e— - l)/r. 



where t^ 2 = X^igCi 2 V 7 "^- ^ or secon d step, we have 
adapted the local detail balance constraint, eq. ([1]), and 
its consequence cq. (|2|) keeping in mind that the corre- 
sponding ratio for the C2 cycles involves only the output 
transition. In fact, the presence of these cycle necessarily 
decreases the power. Moreover, stall conditions (P ou t = 
0) are now reached at a maximum value w ut strictly 
smaller than w m . For the input power, one obtains sim- 
ilarly P in = u; in [(l - e «w-™ in )/ r + + (1 _ e -^)/r 3 +], 
where = X^ec 3 These expressions show that 

both, efficiency r/ = P ou t/Pin and EMP ry* ut , are even 
less universal than for strongly coupled machines. For 
the latter quantity, this fact becomes obvious by look- 
ing at the linear response regime. Expanding the above 
expressions for small < w ou t < w- m <C 1, one finds 



^out 



1/2 



1 + 2(r+/r+ + t+/t 3 + + t^/t+t?) 



+ 0{w m ), 



with all cycle times evaluated in equilibrium. Thus, ?7* ut 
is strictly smaller than 1/2. Still, one should not conclude 
from this result that n* ut necessarily remains bounded by 
1/2 with increasing uij n . 

Concluding perspective.- For soft machines working 
under isothermal steady state conditions, we have in- 
vestigated EMP in the appropriate parameter space by 
varying output power, input power and both. For the 



first two cases, EMP is given by a one parameter func- 
tion whereas in the third case machines universally fall 
into three classes. Our results hold for both unicyclic 
machines and strongly coupled multicyclic ones. Weakly 
coupled multicyclic machines are less efficient and their 
EMP is less universal. Our analytical results and the 
suggested classification provide not only a transparent 
theoretical framework but, in a longer term perspective, 
should also be helpful for designing efficient machines. 
From a theoretical point of view, as a next step, machines 
operating under periodically modulated fields should be 
analyzed similarly for EMP since the extant artificial ma- 
chines typically still require such conditions. 
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